Derivative-free restricted maximum likelihood estimation in animal models with a sparse matrix solver.

نویسندگان

  • K G Boldman
  • L D Van Vleck
چکیده

Estimation of (co)variance components by derivative-free REML requires repeated evaluation of the log-likelihood function of the data. Gaussian elimination of the augmented mixed model coefficient matrix is often used to evaluate the likelihood function, but it can be costly for animal models with large coefficient matrices. This study investigated the use of a direct sparse matrix solver to obtain the log-likelihood function. The sparse matrix package SPARSPAK was used to reorder the mixed model equations once and then repeatedly to solve the equations by Cholesky factorization to generate the terms required to calculate the likelihood. The animal model used for comparison contained 19 fixed levels, 470 maternal permanent environmental effects, and 1586 direct and 1586 maternal genetic effects, resulting in a coefficient matrix of order 3661 with .3% nonzero elements after including numerator relationships. Compared with estimation via Gaussian elimination of the unordered system, utilization of SPARSPAK required 605 and 240 times less central processing unit time on mainframes and personal computers, respectively. The SPARSPAK package also required less memory and provided solutions for all effects in the model.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Derivative-intense Restricted Maximum Likelihood Estimation of Covariance Components for Animal Models

INTRODUCTION Estimation of covariance components by Restricted Maximum Likelihood (REML) fitting an animal model is widely carried out using a derivative-free (DF) algorithm to locate the maximum of the likelihood function. While such algorithms have proven robust and easy to use, they are generally slow to converge, often requiring many likelihood evaluations, in particular for multivariate an...

متن کامل

Running head : REML USING DERIVATIVES Restricted Maximum Likelihood Estimation for Animal Models Using Derivatives of the Likelihood

Restricted Maximum Likelihood estimation using first and second derivatives of the likelihood is described. It relies on the calculation of derivatives without the need for large matrix inversion using an automatic differentiation procedure. In essence, this is an extension of the Cholesky factorisation of a matrix. A reparameterisation is used to transform the constrained optimisation problem ...

متن کامل

The Development of Maximum Likelihood Estimation Approaches for Adaptive Estimation of Free Speed and Critical Density in Vehicle Freeways

The performance of many traffic control strategies depends on how much the traffic flow models have been accurately calibrated. One of the most applicable traffic flow model in traffic control and management is LWR or METANET model. Practically, key parameters in LWR model, including free flow speed and critical density, are parameterized using flow and speed measurements gathered by inductive ...

متن کامل

Approximate Information and Accelerating for High-throughput Heterogeneous Data Analysis with Linear Mixed Models

Linear mixed models are frequently used for analysing heterogeneous data in a 1 broad range of applications. The restricted maximum likelihood method is often preferred to 2 estimate co-variance parameters in such models due to its unbiased estimation of the underlying 3 variance parameters. The restricted log-likelihood function involves log determinants of a 4 complicated co-variance matrix. ...

متن کامل

Lack of Association between Somatotropin Receptor Gene Polymorphism and Birth Weight of Iranian Indigenous Sistani Cattle

The objective of the present study was to determine polymorphism within the promoter region of somatotropin receptor genes in indigenous Sistani cattle (Bos indicus) and associations between this polymorphism and breeding value of birth weight. The pedigree structure was included by considering 1173 animals with 600 progeny birth weight data obtained from a Zhark breeding station in Sistan and ...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Journal of dairy science

دوره 74 12  شماره 

صفحات  -

تاریخ انتشار 1991